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Abstract. We study a random field obtained by counting the number 
of balls containing each point, when overlapping balls are thrown at 
random according to a Poisson random measure. We are particularly 
interested in the local asymptotical self-similarity (lass) properties of 
the field, as well as the action of X-ray transforms. We discover two 
different lass properties when considering the asymptotic either in law 
or on the second order moment and prove a relationship between the 
lass behavior of the field and the lass behavior of its X-ray transform. 
We also describe a microscopic process which leads to a multifractional 
behavior. These results can be used to model and analyze porous media, 
images or connection networks. 

1. Introduction 

The purpose of this paper is the study of a random field obtained by 
throwing overlapping balls. Such a field is particularly well-adapted for 
modeling 3D porous or heterogeneous media. In fact we consider a collection 
of balls in M 3 , whose centers and radii are chosen at random according to a 
Poisson random measure on R 3 x M + . Equivalently, we consider a germ-grain 
model where the germs are Poisson distributed and the grains are balls of 
random radius. 

The field under study is the mass density defined as the number of balls 
containing each point: the more one point is covered by balls, the higher 
is the mass density at this point. From a mathematical point of view, the 
dimension three does not yield any specific behavior, so the study will be 
carried out in dimension d {d > 1). Moreover for d = 2, the number of 
balls covering each point defines the discretized gray level of each pixel in 
a black and white picture. A one dimensional {d = 1) germ-grain model is 
also relevant for modeling communication networks: the germs stand for the 
starting time of the individual ON periods (calls) and the grains stand for 
the 'half-ball' intervals of duration. The obtained process is a counter which 
delivers at each time the number of active connections in the network. 

We have in mind a microscopic model which yields to self-similar macro- 
scopic properties. In order to get this scaling behavior, we introduce some 
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power law behavior in the radius distribution and consider Poisson random 
measures on M. d x M. + with intensity of the following type 



for some 0, which may depend on the location £. The origin of the 'mi- 
croballs' described in this paper can be found in the 'micropulses' introduced 
by Ciosek-Georges and Mandelbrot [Sj with a fixed power 6 in the intensity 
measure. The idea is not new and appeared eighty years ago when S.D. 
Wicksell |19j introduced a first model, the famous 'corpuscles', made of ran- 
dom 3D spheres defined as above. The aim of his study was a stereological 
question. Since then, this kind of model has been extensively deepened and 
extended. We address to ^B] or for many examples of random models 
based on Poisson point process and germ-grain models. Let us also mention 
two recent papers dealing with similar questions. A one dimensional germ- 
grain model with locations (arrival times) uniformly distributed on the time 
axis and intervals lengths (call durations) given by a power law is considered 
by Cohen and Taqqu in [7j. A mixed moving average is performed that sums 
the height of connections and the so-called Poissonized Telecom process is 
obtained. Also similar is the model recently studied by Kaj et al. jlUj : the 
germs are uniformly chosen at random in M. d and the grains are obtained by 
random dilation of a fixed bounded set. In contrast to the quoted models, 
let us point out that inhomogeneity is allowed in our model by choosing a 
non stationary intensity measure © with a non constant power 9 = 6(£). 

This paper is not only concerned with the presentation of a model for 
random media. We also investigate the analysis of the random media mass 
intensity following two ways. On one hand, self-similarity properties are ex- 
plored. More precisely, we focus on a parameter that is supposed to contain 
tangible information on the structure of the media, the local asymptotical 
self-similar index, lass index in short. On the other hand, the action of an 
X-ray transform on the field is explored. This transform is the mathemati- 
cal interpretation for a radiographic process. These techniques are inspired 
from the ones created for Gaussian fields and are still valuable in the Poisson 
context. More specifically we turn to [H] where anisotropic Gaussian fields 
are analyzed by performing X-ray transforms and evaluating lass indices. 
The fundamental aim of these methods is to make a 3D parameter directly 
tractable from X-ray images of the media. 

The notion of local asymptotic self-similarity has been introduced in [2] 
in a Gaussian context and extended to the non-Gaussian realm in |12j . or 
j3j where a general presentation is performed for fields with stationary in- 
crements. The lass index can also be related to other parameters of interest 
as roughness index [2] or Hausdorff dimension pQ. In the area of network 
modeling as well, the notion of self-similarity, at small or large scales, is fun- 
damental and it is highly connected to long-range dependence. The usual 
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self-similarity property requires a scale invariance valid for all scales. This 
is quite restrictive and we will deal with self-similarity properties that are 
fulfilled 'at small scales' only. We introduce a light refinement of the lass 
property of 

Definition 1.1. Let X = {X(x);x € M. d } be a random field and xq £ M. d . 
We call the distribution lass index of X at point xq (fdd-lass index for short) 
the following supremum 

„ fY , J A X0 X(\.)-B(A X0 X(X.)) fdd n n+ \ 
Hf dd (X, x ) = sup <^ a ; — — > as A -> + 

where A Xo X denotes the field of increments at xq : 
A X0 X(x) = X(x + x)- X(x ) 

and — ► means the convergence of the finite dimensional distributions. 

Moreover, when H = Hfdd(X, xq) is finite and when the finite dimensional 
distributions of the centered and renormalized increments \~ H [A XQ X(\.) — 
E(A !ro X(A.))] converge to the finite dimensional distributions of a non- 
vanishing field as A — > + , the limit field is called tangent field at point 
xq (see [Hj). 

When one deals with real data, it is almost impossible to see whether such 
a limit exists in distribution. Therefore we introduce another asymptotic 
self-similarity property, which only uses the second order moment. 

Definition 1.2. Let X = {X(x);x <G M. d } be a random field and xq G M. d . 
We call the covariance lass index of X at point xq (cov-lass index for short) 
the following supremum 

H cov (X, xo) = sup L ; Cov ( A '°* (Aa0 , A ^ Xx 'A _> + V *, 

By analogy with the fdd situation, when the covariance function of X~ H A Xo X(X.) 
for H = Li cov converges to a non vanishing covariance function as A — > + , 
the limit covariance will be called the tangent covariance at point xq. 

Note that the above self-similarity indices are equal for Gaussian fields 
but not in a general setting. Also note that the existence of a tangent 
covariance does not imply existence of the tangent field, and neither the 
converse implication. Actually if H cov is the cov-lass index for X at point 
xo, then for all H < H cov the covariance function of \~ H A Xo X{\.) converges 
to as A — > + . Thus, the finite dimensional distributions of its centered 
version converges also to as A — ► + , and the fdd-lass index Hf^d f° r X at 
point xq - if exists - satisfy LLjdd > Hcov 

Our main results can be summarized as follows: 
- the proposed models provide microscopic descriptions of macroscopic asymp- 
totical self-similar fields which look like fractional -or multifractional- Brow- 
nian motions, depending on the involved intensity measure; 
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- in contrast to the Gaussian case, the covariance lass index and distribution 
lass index are not equal: the first one can be finite whereas the second one 
is infinite, or they can both be finite but with different values; 

- we present explicit formulas that link the lass indices of a field and the lass 
indices of the X-ray transform of it. In particular, when inhomogeneity or 
anisotropy is introduced in the model, it can be recovered through the lass 
indices. 

The paper is organized as follows: the microball model, i.e. the field that 
counts the number of balls covering each point, is introduced in Section 
2. The intensities of the Poisson random measures we will use are spec- 
ified in Section 12.21 A constant power 9 in Q will yield to a fractional 
microball model which is stationary and isotropic. A non-constant power 
will yield the multifractional model. We also introduce in Section 12.31 the 
X-ray transform. Section 3 is devoted to the self-similarity properties of the 
microball model and its X-ray transform. Theorem 13. II deals with the frac- 
tional model and Theorems 13 .41 and 13. 61 deal with the multifractional model, 
where = is a smooth, or a singular function respectively. We also 
compare our results to homogenization results in Section f3.3l The proofs of 
Theorems 13.41 and 13.61 are detailed at the end of the paper. The covariance 
lass properties are proved in Section [51 and the distribution lass properties 
in Section [SJ 

2. The microball model and its X-ray transforms 

2.1. The random grain model. The model is built by considering the 
superposition of balls B(£, r), where £ is a point in W d and r > is the radius 
of the ball. As in we want to study the mass distribution generated by 
a family of balls B(^j,rj), with random location £j and random radius r,-. 
We assume that (£j,Fj) are given by a Poisson point process with intensity 
v(d£,dr), where v is a non-negative cr-finite measure on R d x R + . 

For each x € M rf we are naturally interested in the number of balls B(£, r) 
that contain the point x, given by 



Moreover, in that case, we can represent the field through a stochastic inte- 
gral with respect to a Poisson measure N with intensity v, as 



Let us also describe the intuitive scenario we have in mind for the one dimen- 
sional case (d = 1). We look at the 'half-ball'-interval + r) as the ON 
period of a single call and the number of connected users at time x is equal 



#{j;x e B(£ j} rj)} = l B ^ r ^(x). 



j 



Such a field is well defined as soon as 
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to the one dimensional integral -if exists- / KxK + l[££+ r )(x)N(d£, dr). Since 
the integrand behaves like l B ^ r ^, we will not distinguish anymore 

the special case d = 1. 

2.2. The microball model. As in j^j, we assume that the radii of such a 
random grain model follow a power law in (0, 1). The exponent of the power 
law can either be constant (fractional case) or can depend on the center of the 
ball £ (multifractional case). More precisely we consider intensity measures 
on R d x R+ with special shapes v m and Uh described below. 

For m > we define the fractional intensity measure v m on R d x R + as 

(3) v m (dt, dr) = r- d - 1+2m l (0A) (r)d^dr . 

For a function h on R rf such that essinf h > 0, we define the multifractional 
intensity measure Vh on R rf x R + as 

(4) Md^dr) = r- d - 1+2h ^l {0A) (r)d(dr . 

It is straightforward to see that these kind of measures satisfy (|2*|). This 
allows the following definition. 

Definition 2.1. Let h be a function on M. d such that essinf h > 0, and let 
Nh be a Poisson random measure with intensity . The field defined on M. d 
as 

X(x) = [ l m>r) (x)N h (dt,dr) 

is called a microball model with index h. 

Note that the microball model has moments of all order. In particular, 
its mean value is given, for each x € R d , by 

B(X(x)) = [ l B{M (x)v h ((%,dr). 

Moreover, by the isometry of the Poisson measure, its covariance function 
is equal to 

Cov(X(x),X(x')) = / l B K r )(x)l B ^ r )(x')i> h (d£,dr), 

JR d xR+ 

for all x, x' G R d . 

The self-similarity properties of the microball model that we will study 
in the sequel deal with the local behavior of X. We now compute the 
increments of X and analyse their moments in the following lemma. In 
what follows we write for xq,x £M. d 

(5) A X0 X(x) = X(x + x) - X{x Q ) = / i/,(x,Z-x ,r)N h (d£,dT), 

JR d xR+ 

where 

(6) il>(x,Z,r) = l B ^ r) {x) - lB(f,r)(0) = 1 \"-i\<r<\& - ^-ISI^k-CI' 
and l.l denotes the usual Euclidean norm. 
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Lemma 2.2. Let m G (0,1/2). There exists a constant C(m) G (0, +oo) 
such that, for all p G (0, +oo) and x E M, d , 

|V(x,£,r)| p r- d - 1+2m drdi = C{m)\x\ 2m . 

'xl+ 

Proof. Note that |^(a;,^,r)| p = \ip(x,£,r)\ = l\ x -£\< r <\t\ + 1 [el <r-< ^| • 
Hence, for all x, £ G 

^(x^,r)P r- rf - 1+2m dr = ^A_i |€|<|ie _ €| (|£r d+2m _ \ x _ ^-d + 2m^ 

and for < m < 1/2, the function £ i-> |£|- d + 2m — \ x — £\~ d + 2rn j s integrable 
on R . The conclusion is obtained by rotation invariance and homogeneity. 

□ 

2.3. X-Ray transform. One motivation for this paper is to describe, model 
and analyze heterogeneous media. We have in mind the possibility to esti- 
mate a macroscopic 3D parameter through X-ray images. By this method, 
it will be possible to get an analysis of the media without entering the me- 
dia (non-invasive method). In this section the mathematical tool associated 
with X-ray images is presented and tested on the microball model. We 
assume that d>2. 

Following the usual notation (see JS] for instance), for a direction a E 
S^ 1 = {x G M. d ; \x\ = l}, the X-ray transform in the direction a of any 
function / G L 1 (M d ) is given by 



y G< a > ± / f(y+pa)dp 
Jr 



where <a> i := {ieR li ;i'« = 0}, and ■ denotes the usual scalar product 
on M. d . We want to define, in the same way, the X-ray transform of a 
microball model X. Unfortunately, the realizations x G M d i— > X(x,u>) 
do not belong to L 1 (R d ). We will therefore work with the windowed X- 
ray transform defined through a fixed window p. We assume that p is a 
continuous function on M with fast decay, namely 

(7) V7V G N, 3CV, Vp G M, \p(p)\ < C N (1 + \p\)~ N . 

For any function / G Lj oc (M. d ) with slow growth, we define the windowed 
X-ray transform of f in the direction a to be the map 



(8) y G< a > 



Pafiy) ■= / f(y + pa)p(p)dp. 

JR 



Our aim is to apply such a transformation to a microball model. We use 
the properties of the Poisson random measure and the fact that one can 
define the field 

\ I ( / 1 B(t,r)(y + P a )p(p)dp)N h (d£,dr)) y G< a > x \, 

lJR d xR+ \JR / J 
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as soon as the integrand 

'Pa (lB(£,r)(0) (y) = / W.r) {V + P&)p(p)dp 

Jr 

is integrable with respect to ^^(d^dr) for each y G< a > . But this is 
straightforward using © and (J7J). This allows us to state the following 
definition. 

Definition 2.3. Let h be a function on M. d such that essinf h > and let 
Nh be a Poisson random measure with intensity v^. Let X be the microball 
model with index h, X(-) = J Rt e xK + ^-B(^r)(,')^h(dS,, dr). 
Let p be a continuous window that satisfies Q). For a £ S d ~ l , we call 
the windowed X-ray transform of X in the direction a, the field defined on 
< a > ± by 

V a X(y) = [ V a (lB«,r)(0) (y)N h (d£,dr). 

jR d xR+ 

Note that for (£, r) in R d x M + , Cauchy-Schwarz inequality leads to 

\V a {l B ^ r) (-))(y)\ 2 < 2(r 2 -\y-7r a± (0\ 2 )T [ l B ( € ,r) (y + pa) p 2 (p)dp 

Jr 

< 2r / W, r ) (V + P«) P 2 (p)dp , 
JR 

where 7r Q ± denotes the orthogonal projection on < a > . Thus, for all y G< 
a > x , the function "P a (l B /^ r \(.)) (y) belongs to L 2 (R d x R + , i>h{di, dr) and 
admits a second order moment. 



As for the microball model itself, we are particularly interested in the local 
behavior of the X-ray transform and we need to estimate its increments. 

For yo, y G< a > , using (JSJ) we get 
(9) 

A yo V a X(y)=V a X(yo + y)-V a X(y )= [ G p (y,£ - y ,r)N h {d^,dr) 

JR d xR+ 

where 

(10) G p (y,€,r)= / i/)(y,Z - pa,r)p(p)dp . 

Jr 

Note that the above integral is well defined for any bounded function p. 
In the special case where p = 1, we write G instead of G\ and for j/,7 in 
< a >^ and r in M + , a simple computation gives 

(11) G(y, 7, r) = G 1 (y, 7, r) = (r 2 - |y - 7 | 2 )f - (r 2 - | 7 | 2 )f, 

where, as usual, t + : = max(0,t) for all t 6 R. 

The next lemma provides upper-bounds for the integral of G(y, .). 
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Lemma 2.4. For m € (0, 1/2) and y €< a > , 

G{y,.) € L 2 (< a > x xl+ r- d - 1+2m d 7 dr) . 
Proof. For m 6 (0, 1/2), y G< a >-*- and A > 0, on one hand 



G{y,l,rfd-( < r 2 



<a> A 



<a> J 



l \y~7\<r 



+ l| 7 |<r) ^1 



(12) < C r d+1 . 

On the other hand, a change of variable gives, for y ^ and r > 0, 
■[„|N -<*-!> 



G(y,j,r) dj 



<a> J 



1 



\y\ 



<a>- 
2\ V2 



2r 



7 + 



1/2 



d^ 



The next lemma provides an upper bound for the last quantity, which leads 
to 

r 2r 
(13) / G(y, 7, r) 2 d 7 < C\y\ 2 r d ' 1 ln(2 + n ) . 

J< a >^ \y\ 

Since m € (0, 1/2), inequalities (fT2*|) and (|TT?|) conclude for the proof. □ 

Lemma 2.5. Let n £ N*. There exists a constant C > suc/i i/ia£ /or a// 
direction e E <S n_1 and a// r > 0, 

2 



/ 2 1 1 2 \ 1/2 / 2 1 1 
(r — \x — e\ ) / — (r — |x + e 



2\l/2 



< Cr n ln(2 + r). 



Proof. For n = 1, we have to prove that there exists a constant C such that, 
for r > 0, 

+1 (r 2 - (x - l) 2 ) 1 ' 2 - (r 2 - (x + 1) 2 )+ /2 2 dx<Cr ln(r + 2). 



This is an easy consequence of the fact that the function that we integrate 
is bounded by 4r for x G [r — 1, r + 1], and by 16r 2 ((r — l)(r — x + 1)) _1 for 
x € [0, r — 1] when r > 1. 

In the general case (n > 1) we write x = x' + x"e with x' = vr <e> ±(x) and 
1" £ I. From the one-dimensional case, for x' £< e >^ , 

2 

dx" 



2 M|2_| x // + 1 |2 ) l/2 



( r 2_| x '|2_| x //_ 1 |2 ) l/2_ (r 

< C{r 2 -\x'\ 2 ) 1 l 2 ln(r + 2). 



But 



<e> J 



JO 



(r 2 -|x'| 2 )^ dx' = r r 
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Finally, we can change the constant C such that 

[r 2 -\x- el 2 ) 1 / 2 - (r 2 -\x + e\ 2 )]/ 2 dx < Cr n In (2 + r). 

□ 

3. Lass properties 

The section is devoted to the study of self-similarity properties of the 
microball models and their X-ray transforms. 

3.1. The fractional microball model. Let d > 1. Let m > and recall 
that the fractional intensity measure v m on M. d x R + is given by 

v m (dt,dr) = r- d ' 1+2m l m (r)didr 

and the fractional microball model by 

X = \ [ l B ^ r) (x)N m (d^,dr); x G R d ) , 

lJR d xR+ ) 

where N m is a random Poisson measure of intensity 

First, let us remark that the choice of v m as intensity means that the 
center and the radius are thrown independently. Hence, since the centers 
are uniformly distributed on the state space M. d , the fractional microball 
model X is isotropic and stationary, i.e. for every rotation R centered at 
in R d and for all x G R d , 

X o R f = X and X{x + .) f = X . 

In the following, we will assume that m < 1/2 so that Lemma 12 . 21 applies. 
Then, using (JSJ) and the stationarity of X we get for x, x' G 

E (X(x) - X{x')f < C{m)\x - x'\ 2m . 

Thus, the field X is mean square continuous. 

Furthermore, by the correlation theory of stationary random fields (see 
[ZDj for example), there exists a finite positive Radon measure a such that 

Cov(X(x),X(0)) = / e- ix ^d(r(0- 



By computing the inverse Fourier transform we obtain that a is absolutely 
continuous with respect to the Lebesgue's measure. Moreover, the spectral 
density of the fractional microball model X is given by 

'd/2\ 



where is the Bessel function (see p. 406 for example). Let us recall 
that 

d/2 

Jd ^ = 2d/ 2 T{1 + d/2) (i + o(M 2 )),^o, 



10 HERMINE BIERME AND ANNE ESTRADE 

and 

Jd/a W = ]ff s (cos (s - + o (s" 1 )) , +00. 

Hence, the previous spectral density has a power law behavior and lass prop- 
erties for the associated field are expected. This point of view is detailed in 
01 section 2.2.2. 

We first investigate the covariance lass property of X. Recall that X 
is stationary and therefore the behavior of the increments of X around any 
xq G M. d does not depend on xq. We have to deal, for x, x' G M. d and A — > + , 
with 

Cov(A X0 X(\x), A X0 X(Xx')) = [ t/;(\x, £, r)^(Ax', e, r) r - d - 1+2m drd( , 

JM d x(0,l) 

where ip is given by (jHJ) • By homogeneity and integrability of if) (see Lemma 
I2.2|) it is easy to establish that, for all x, x' € R d , 

lim A ^ 0+ \- 2m [ ^(\x,tr)^(\x',£,r) r- d - 1+2m drd£ 

(14) = ( 7p(x,t,r)7p(x',£,r) r- d - 1+2m drd£ . 

The cov-lass property of X follows from 1)14|) . 

Now, let us assume that d > 2 and choose q S 5 . We look for the 
cov-lass property of V a X, the windowed X-ray transform of X. Let us recall 
© and (flT)]) : for any yo and y in < a > , 

A yo P Q X(y)= / G p {y,ti-y Q ,r)N m (d£,dr) 

where 

G p (y,€,r)= / - pa,r)p(p)dp . 

Jr 

It is straightforward to see that V a X is also stationary so that the cov-lass 
properties will not depend on the point yo- 

Writing any £ G M. d as £ = 7 + ta with 7 G< a >^ and t £ R and 
performing a change of variable (translation-dilation) in the integral that 
defines G p , we obtain 

(15) Cov(A yo V a X(Xy), A yo V a X(Xy')) 

= A 1+2m / G p(t+A .) (y, 7, r)G p(t+A0 (y', 7, r) r^^ 2 ™ d 7 dtdr, 

i<a> 1 xlx(0,A- 1 ) 

with p(t + A.) denoting the window p 1— > p(t + Ap). Note that for r € R + , 
teM and y, 7 G< a > ± , 

Gp(t+A.)(y,7^) ^ G(y,-y,r)p(t) . 
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Assumption @ on/) and Lemma 12.41 allow to conclude that 
lim X~ 1 - 2m Cov(A yo V a X(Xy),A yo V a X(Xy')) 

A — *0~^~ 

p(t) 2 dt] [ G(y, 7, r)G(y', 7, r )r- d - 1+2m d 7 dr. 

: / J<o> 1 xt+ 

We now state a theorem including the notions of local asymptotical self- 
similarity index given in Definitions II. II and II. 21 The fdd-lass properties are 
proved in Section as a consequence of Theorem 13.41 in the special case 
where the multifractional index h is constant equal to m. 

Theorem 3.1. Let d > 2. Let X be the fractional microball model with 
index m S (0,1/2) and let V a X be its X-ray transform in the direction 
a e S d -\ 

• At any point xq £ M, d , the cov-lass index of X is equal to m and 
the fdd-lass index of X is equal to +00. Moreover the covariance 
of X~ m A Xo X(X.) converges, up to a multiplicative constant, to the 
covariance of a fractional Brownian motion of index m. 

• At any point yo £< a > , the cov-lass index and the fdd-lass index of 
V a X are equal to m+1/2. Moreover the covariance and the finite di- 
mensional distributions of A"™- 1 / 2 (A yo (7> a X)(\.) - E (A yo (V a X)(X.))) 
converge, up to a multiplicative constant, to the corresponding ones 

of a fractional Brownian motion of index m + 1/2. 

Remark 3.2. The first point of Theorem \3.1\ is still true in the one- dimensional 
case (d = 1). 

Let us comment this result. 

In one dimension, this result describes the small scale behavior of the 
number of active connections in a communication network: the covariance 
is locally asymptotically self similar and behaves like a fractional Brown- 
ian motion covariance. More generally, the same is observed in the multi- 
dimensional case. Hence the fractional microball model provides a micro- 
scopic description of a random media which behaves, up to the second-order 
moment, like a fractional Brownian motion. 

The second point of Theorem 13. H is very interesting from a practical point 
of view. The one-to-one correspondence between the lass index of X and 
the lass index of V a X allows the estimation of the 3D lass index through the 
analysis of the media radiographic images. The fractional microball model 
is thus relevant to model isotropic and stationary media. 

Following a widespread idea ( ^1] , [3] , [H] ) > we now consider the multifrac- 
tional microball model, where m is replaced by a function that depends on 
the ball location. 



3.2. The multifractional microball model. In this section we study the 
lass properties of the multifractional microball model, associated with the 
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multifractional intensity measure u^, given on M. d x IR + by 

v h (di,dr) = r- d - l+2h ^l {0 , 1} (r)dtdr , 
where h is a function on M. such that < h(£) < 1/2 . 

First, let us remark that the multifractional microball model is not sta- 
tionary nor isotropic when h is not constant. 

Let us deal with the cov-lass properties of the multifractional microball. 
One has to study, for all xq, x, x' G R , and A — > + , the asymptotic behavior 
of Cov (A X0 X(Xx), A XQ X(Xx')). By a change of variables, with tp given by 
®, 

(16) Cov (A X0 X(Xx), A XQ X(Xx')) 

= [ A 2,l(:ro+A « ) V'(2;,e,?')V'(^,C^)r" d " 1+2 ' l(:i ' 0+A0 ^r. 

Jr^x^.A" 1 ) 

Similarly, when d>2 and a G S^ 1 , for all yo,y,y' G< a > J -, Q15|) becomes 



(17) Cov(A yo V a X(Xy), A yo V a X(Xy')) 

= A / G p{t+x ^(y,-f,r)G p{t+x ^(y',-f,r) 

J<a> J -xRx(0,\- 1 ) 

x (XrfHyo+ta+x^^d-i djdtdr 

In order to get cov-lass properties for both X and its windowed X-ray 
transform, further assumptions on h have to be made. We are mainly inter- 
ested in two kinds of functions. The first kind deals with the case where h 
is smooth on R . It is linked with the multifractional Brownian motion 14 , 
0| , obtained by substituting the Hurst parameter H by a Lipschitz function 
on the state space. The second class of function is when h is singular at the 
point 0, defined on M. \ {0} by homogeneity, ie h(X£) = h(£) for all A G M*. 
This follows the point of view of |Hj to get anisotropic generalizations of the 
fractional Brownian motion. 



Let us recall the definition of a /3-Lipschitz function. 

Definition 3.3. Let (M, dja) be a metric space and [5 G (0, 1]. A function 
f : M — > R is called (3-Lipschitz on M if there exists C > such that 

Vx,y G M, d M {x,y) < 1 \f(x) - f(y)\ < Cd M (x,yf. 

3.2.1. The smooth case: let us first study the case of a /3-Lipschitz func- 
tion on M. d . By continuity of h around xq G R d , we get from (|16|) that 
H cov (X, xq) = h(xo), using Lemma l2~2l A rigourous proof of this statement 
is given in Section [SJ Moreover, the continuity of h around yo + tot, for 
yo G< a >- L and t G R, applied in (jTTj) and Lemma EU will imply that 
H cov (V a X,y ) = m(a,y ) + 1/2, where 

(18) m(a, yo) := inf h(yo + ta). 
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These observations yield to the next theorem. A detailed proof is given in 
Sections El and El We denote meas the Lebesgue's measure. 

Theorem 3.4. Let d > 2. Let h be a (3-Lipschitz function on M d such that 
< h < 1/2. Let X be the multifractional microball model with index h and 
let V a X be its X-ray transform in the direction a G S . 

• At any point xq G M. d , the cov-lass index of X is equal to h(xo) and 
the fdd-lass index of X is greater than min(l, f3 + 2/i(xo)). Moreover 
the covariance of \~ h ( x °} A X0 X(X.) converges, up to a multiplicative 
constant, to the covariance of a fractional Brownian motion of index 
h(x ). 

• At any point yo G< a >~ L , the cov-lass index and the fdd-lass in- 
dex ofP a X are equal to m(a,yo) + 1/2. Moreover, when meas({t G 
R*! h(yo+ta) = m(a,yo)}) > 0, the covariance and the finite dimen- 
sional distributions f\- m ^) +l l 2 (A yo (V a X)(\.) - E (A yo (P a X)(\.))) 
converge, up to a multiplicative constant, to the corresponding ones 

of a fractional Brownian motion of index m(a, yo) + 1/2. 

Remark 3.5. The first point of Theorem \3.J\ is still true in the one- dimensional 
case (d = 1). 

We remark that the cov-lass index of X at point xq G M rf is equal to h(xo), 
as for the multifractional Brownian motion with index h, and this justifies 
the name multifractional. Also let us just point out that these results apply 
to the fractional case, where h is constant. In this case we have shown that 
the fdd-lass index for the microball is equal to +oo. This is a consequence of 
the stationarity of the model. The first order moment of the increments is 
then equal to and this is no lomger the case when h is not constant. From 
a certain point of view, for the microball model, the fdd-lass index deals 
with the regularity of the first order moment and the cov-lass index with 
those of the second order moment. Finally, compared to the fractional case, 
we still have an additive factor of 1/2 for the lass indices of the windowed 
X-ray transforms, but here only the infimum of h along straight lines can 
be recovered. 



3.2.2. The singular case: let us now consider h to be an even, /3-Lipschitz 
function on the sphere, defined on M. d \ {0} by homogeneity. Of course, 
this case is only relevant when d > 2 because otherwise it turns to be the 
fractional case, with h constant. The singularity of h at point makes this 
point a very special one. The balls are thrown from and their numbers 
and sizes only depend on the direction along which they are thrown. To 
distinguish the microball model associated with such a singular h from the 
multifractional one, we will call it the star microball model. Let us remark 
that the /3-Lipschitz assumption on h means that there exists C > such 
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that for all x G M d \ {0} and f G M d , when j£| < 1 
(19) |M^o + -h(xo)\ <C\x \-P\Hf. 

Then, the lass properties of the star microball model, respectively its 
X-ray transform in the direction a G S' rf_1 , at any point xq G M. d \ {0}, 
respectively yo G< a >-*- \{0}, are the same as for the multifractional mi- 
croball model, respectively its X-ray transform in the direction a, given in 
Theorem 13.41 Thus, the next theorem will only deal with the lass properties 
around 0. Let us remark that, in that case, by homogeneity of h, the expo- 
nent in ()16|) is equal to \ h ^' for all £ ^ 0. Then, denoting by m G (0, 1/2) 
the minimum of h on S l , by Lemma |2.21 we will get H cov (X, 0) = m. 
Moreover, in (fT7|) . since /i is continuous around to and h(ta) = h(a), for 
t G E*, using Lemma l2~H we prove that H cov (V a X,0) = h(a) + 1/2. 
Let us state the different lass properties of the star microball model. 

Theorem 3.6. Let h be an even (3-Lipschitz, non constant, function on 
S^ 1 such that < h < 1/2. Let X be the star microball model with index h 
and let V a X be its X-ray transform in the direction a G S"^ -1 . We consider 
lass properties at point 0. 

• The cov-lass index of X is equal to m and the fdd-lass index of X is 
equal to 2m. Moreover, when meas({h = m}) > 0, the covariance of 
\~ m AoX(\.) converges to 7 m with 

j m (x, x>) = ( l m=m ^{x, e, r)^{x',i, r)r- d - 1+2m d(dr , 

JR d xR+ 

for x, x' G M. d , while the finite dimensional distributions of 

\- 2m (A X(A.) - E (A X(A.))) converge to the deterministic field 

Z m , with 

Z m = \~\x\ 2m [ l h(€ )=m^(ri -^,r)r- d - 1+2m d^dr ;x£l4, 
I jR d xR+ \ x \ ) 

where ip is given by JSJ). 

• The cov-lass index and the fdd-lass index ofV a X are equal to h(a) + 
1/2. Moreover, the covariance and the finite dimensional distribu- 
tions of X~ h ^ +1 / 2 (A (V a X)(X.) -B(A (V a X){\.))) converge, up 
to a multiplicative constant, to the corresponding ones of a fractional 
Brownian motion of index h(a) + 1/2. 

The proof is given in Sections and El Let us remark that for the star 
microball model, there exists both a cov-lass index and a fdd-lass index and 
that the later equals the double of the former. This multiplicative factor is 
typical for the Poisson structure proved by the following exercise. 

Remark 3.7. Let (X n ) be a sequence of Poisson random variables. Suppose 
there exists some H > and v > such that Var(n H (X n — E{X n ))) tends 
to v when n tends to +oo. Then n 2H (X n — E(X n )) tends in distribution to 
—v. 
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Moreover when {h = m} has positive measure, the tangent field at is 
deterministic and not zero, hence does not have stationary increments. This 
is worth to be noticed and linked to a result of Falconer [B|, which states 
that at almost all points the tangent field -if it exists- must have stationary 
increments. Hence the point appears as an 'exceptional point' (see ^H] for 
other examples of exceptional points). 

Finally, let us point out that the tangent field of the X-ray transform, 
when it exists, is Gaussian, even a fractional Brownian motion, whereas the 
tangent field of the star microball model was deterministic. This justifies, 
from a mathematical point of view, modeling radiographic images by fBm, 
even when the media under study is far from being of this type (see [S] for 
an experimental study). 

3.3. Comparison with homogenization results. There are different ways 
to consider self-similarity at small scales, depending on which part of the 
signal the scaling acts. Instead of performing a scaling on the increments 
lag, as it is performed in Sections 13.11 and 13.21 we act on the radius of the 
balls as follows. Suppose we zoom and consider the balls B(£,r/e) instead 
of the balls B(£,r), where the (£, r) are randomly chosen by the Poisson 
random measure Nh, and we let e tend to 0. Denoting by X £ the associated 
field 

X £ (x) = f l mr/e) (x)N h (d£,dr) , x e R d , 

JR d xR+ 

we look for normalization terms n(xo) such that e n ^ x °\A XH X £ — E(A X() X e ) 
converge in distribution to a non degenerate field. Note that the field X s 
can also be considered as a microball model (see Definition 12. 1J) associated 
with a Poisson measure with intensity 

v £ h (d£,dr) = e-^^r-^+^h^-i^d^dr . 

Actually this procedure is nothing but homogenization and is close to the 
thermodynamical limit investigated in jHj or the scaling limit in jlUj . Similar 
computations as for the previous theorems yield 

• if < m := min/i < 1/2 and the set {£; h(£) = m} has positive measure, 
then the normalization term n(xo) is equal to d/2 — m for all xq G W 1 

• moreover if h is constant equal to m, then the limit field is a fractional 
Brownian motion with index m. 

4. Conclusion 

We propose to model -from a microscopic point of view- the mass intensity 
of a porous media or the number of connected custumers in a network with a 
non-Gaussian field, which presents a macroscopic (multi)fractional behavior. 
The rich structure of Poisson point processes allows us to reach this goal and 
also to perform explicit computations as in the Gaussian case. In order to 
keep the model as intuitive as possible, we did not try to produce more 
general fields. The Poisson structure can obviously be exploited further on 
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by considering more general integrators with respect to the Poisson measure. 
Replacing the indicator function 1b(£,t-) i n Definition 12 . 1 1 bv a more general 
one will, for instance, give the possibility to model granular media with non 
spherical grains. Another model for porous media can also be built up from 
a collection of random spheres which correspond no more to grains, but to 
pores or bubbles. By this way, one will get a {0, l}-valued field and leave 
the linear context. 

5. Proofs of the cov-lass properties 

Let h be a function defined on R d or R d x {0}, with < h < 1/2. 
In this section we will give rigorous proofs for the cov-lass properties of the 
microball models and their windowed X-ray transforms. 

5.1. Cov-lass properties of the microball model. Let xq G M. d . For 
He (0, 1) and x, x' G M d , recall from (HHJ) that 

r?(W):=Cov(^^,^W>) 

= f A- 2(H - h(xo+A€)) V'(x,e,r)V'(x , ,e,r)r- d - 1+2 ' l ^ 0+A ^^(ir. 
JR' i x(0,A- 1 ) 



a. The smooth case: we assume that there exists C > such that, for 

iei<i 

(20) \Hx + 0-h(x )\<C\tf. 
We first establish that 

(21) T h x ixo) {x , x, x') r^°) (x, x') , 
where 



F h ^ x °\x,x') 



R d xR+ 



ip(x, £, r)il>(a/,£, r)r- d - 1+2h( - X0 U£dr 



and then prove that r h(xo) is the covariance function of a fractional Brownian 
motion with Hurst index h(xo). 

Since h(xo) < 1/2, Lemma 12 . 21 indicates that the function ip{x, .) belongs 
to L 2 (R d x M. + ,r~ d ~ l+2h ( xa '> d£,dr) . By Cauchy-Schwarz inequality, in order 



to prove <|2~T|) . it is enough to prove that the difference /(A), given by 
(22) 

/(A) = f ^(x,£,r) 2 (Ar)- 2 (' l ( :Co )- /l ( a; o +A «)) - 1 r - d - 1+2h ^ d£dr , 

tends to as A tends to + . Let us remark that, by Q2UJI . for some positive 
p and q G (0, 1) to be fixed later, when A p+1 < Ar < 1 and Aj£| < A 1-9 < 1, 
one can find C > such that 



(23) 



(Ar) 



-2(h(x )-h(x +\i)) 



1 



< C(A|£|f [ In A[. 
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Without loss of generality, we can assume /3 < 1 — 2h(xo). In this case we 
will prove that, for A small enough compared to x, 

(24) /(A) < CA^mAl, 

with C > 0. We split the integral into 



/ -/ 

itfxfO.A" 1 ) JB 



+ 



S(0,A-'3)x(Ap,A- 1 ) 



/ +/ 

JR d x(0,AP) JB 



x(0,Ap) JB(0,A-?) c x(0,A- 1 ) 



By (|23j) , the same kind of arguments as in the proof of Lemma 12.21 yield to 

< CA^|x|^ +2 ^ D )|lnA|. 



'B(0,A-9)x(Ap,A- 1 ) 

Moreover, for £ € B(0, A~ q ) c or r E (0, A p ), we use the following inequality 



(Ar 



,-2(/i(x )-ft(xo+AO) 



1 



< 2(Ar 



.2(m— /i(xo)) 



which holds since /i(xq + A£) > m and Ar < 1. Then, on one hand, 



! x(0,AP) JR d x(0,\P) 
^ ^r^2(m-fe(3;o))+2nip 

On the other hand, for A < (4(1 + 



■\£\<r 



-d-l+2m 



S(0,A-9)c x (0,A- 1 ) 



-d— l+2m 



< C|x|A 2(m - /l{a;o)) 

^iei>A-8 

< (7| x |_\ 2 ( m -' l ( :r o))+i2(l-2m) 

Finally, it suffices to choose p and g such that 

3 + 2(h(x ) - m) 3 + 2(h(x )-m) 

v > o and 1 — o < g < l- 

2m 1 — 2m 

This concludes for the proof of the convergence of r h x {xo \x ,.) to r^°). it 
remains to show that T h( - X °^ is the covariance of a -up to a constant- fractional 
Brownian motion with Hurst index h(xo). A straightforward computation 
gives 

(il>(x,£,r) - i){x',Z,r)) 2 r -d-Wh(xo) d£dr 



/R d xR+ 

This allows to write T h ( x °^ as 



T h ( xo \x, x') = i (v(x) + w(x') - u(x - x')) 
where, = T h ^(x, x) = c\x\ 2h ^ . 

This proves the covariance part of the first point in Theorem 13 .41 
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b. The singular case: we now assume that h is given by an even (3- 
Lipschitz function on the sphere, defined on M. d \ {0} by homogeneity, and 
deal with the case xq = 0. 

b-i) Assume that {£ £ M. d \ {0}; h(£) = m} has positive measure. 

Thus, using Lemma \2.2l Lebesgue's Theorem gives that T™(0, x, x') tends 

to 

(25) / l m=m ^(x,^r)x/;(x',C,r)r- d - 1+2m dCdr := 7 m (^'), 

JR d xR+ 

which does not vanish by assumption. 

b-ii) Assume now that {£ € M. d \ {0}; /i(£) = m} has measure 0. 

Since j m (x, x') = by (|23|) . the covariance if (0, x, x') tends to for H < m. 

On the other hand, for H = m + 2e with e > and for A in (0, 1), 

if {x, x) > A- / l m <m+e^(x, e, rf r- d ' l+2h ^d^dr , 

JR d x(0,l) 

where {£ G M d \ {0}; h(£) < m + e} has positive measure. So the above 
quantity tends to infinity when A tends to + . Hence the exponent m is 
proved to be the cov-lass index for X at xq = 0. 

The cov-lass properties asserted in the first part of Theorem 13.61 are now 
proved. 

5.2. Cov-lass properties of the windowed X-ray transforms. Let a G 

S" i_1 and yo S< a >-*-. Let us denote by T^(yo, .) the covariance function 
of the field \~ H Ay V a X(\.). Then for y and y' in < a > ± , by (|T7|) 

if iyo,y,y') = \~ 2H+l \ G P (t+\.)(y,i,r)G p{t+x .)(y',-/,r) 

J<a> ± xRx(0,\- 1 ) 

x (\ r f h ^+ ta +y*) r - d - 1 d^dtdr. 
Let us write G p (j+a .) as the integral given by ((TU|) 

rf {yo,y,y') = \~ 2H+l / V"(y,7 -pa,r)p(t + Xp)(y,^,r) 

J <a> ± xRx(0,\- 1 )xR 

x G p(t+A . ) (y / , 7 ,r)(Ar) 2fc ^ +to+ ^)r- d - 1 d 7 cftdrdp. 
Two changes of variables allow to write if (yo,y,y') as 

\~ 2H+l f ^(y, e, r>(t)G p(t+ A.) (y' } £, x (Ar) 2/l ( A « +to+ ^)r- d - 1 d£drdi. 

a. T/ie smooth case: we assume that /i is /?-Lipschitz. By the same 
arguments as in the proof of (|24j) . using the fact that G p < 2 H/oHoo |G| , for A 
small enough compared to y and y' , 

2h(\Z+ta+yo)-2h(ta+y ) _ ^ 




^(y,C^)G , p( t+A .)(y / ^> r ) 
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(26) x r -d-i+2h(t*+yo) d ^ dr < CA /3| ln ( A )| ; 

with C > and (3 < 1 - 2h(ta + y ). Thus we write (\ r fH*>Z+t<*+yo) as 

^ r ^2h(\£+ta+y )-2h(ta+y ) x ^ r ^2h(ta+y )-2m(a,y ) x ^ r ^2m(a,y ) 

where m(a,yo) is given by (|18|) . For almost every i € R, the second factor 
tends to 1-h(yo+ta)=m(a,y )- Since m(a, yo) < 1/2, let us recall that by Lemma 
H3]the function G(y, .) belongs to L 2 (< a > x xR+, r -d-i+2m(a,y ) d7dr y 
We use Lebesgue's Theorem and (|2*6*|). to obtain the following asymptotics: 
rf a > y ° )+1/2 (y ,y,y') tends to 

(jf l%o+*a)=rn(a,«o)P 2 (^) x r m ( Q ^) +1 / 2 (y,y'), 

where 

(27) T H (y,y'):= [ G(y, 7 ,r)G{y' , 7 ,r)r- d - 2+H d 7 dr. 

The identification of as the covariance of a fractional Brownian motion 
defined on < a >-*- with Hurst index H = m(a, yo) + 1/2 is straightforward 
following the same arguments as in the part a. of 15.11 

If {t S R;/i(yo + ta) = m(a,yo)} has positive measure, then we have 
finished with the proof. Otherwise, we proceed in the same way as in part 
b-ii) of 15.11 The cov-lass property of the second part of Theorem 13.41 is 
established. 

b. The singular case: we assume that h is /3-Lipschitz on S^ -1 with 
< (5 < 1 and look at the cov-lass property around 0. By homogeneity, the 
/3-Lipschitz condition is replaced by 

\h(ta + x) - h(a)\ < Ct^lxf. 

Hence the upper bound of (|26j) is now given by 

C (|i|- (1 ~ 2m) + |r (1 - 2/l(Q)) ) \ p \ ln(A)|. 

Lebesgue's theorem still applies to get the result enounced in the second 
part of Theorem 13.61 

6. FDD-LASS PROPERTIES 
We now prove the finite dimensional distribution lass properties. 
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6.1. Fdd-lass properties of the microball model. Let us denote X = 
X — E (X) the centered version of X. For notational sake of simplicity, we 
will only consider the limit in distribution of \~ A X0 X(\x) for a fixed x 
in M. d instead of a random vector (A~ H A Xo X(\xj))i<j<. n . The general case 
follows along the same lines. 

For H > 0, x € R d and t G R, let 

/ A X0 X(Xx)\ 
Eexp I — -i-^ — - I = exp&(H, \,x ,x,t) 

where &(H, A, xq, x, t) is given by 

I «A-^(Ax^-x ,r) _ x _ riA - H V(Ax, e - x , r)) r) . 

! xM+ V 7 

A change of variable yields 

<Z>(H,\,x ,x,t) = [ (e itx ~ H ^ x ^ r) -l-it\- H ^{x,i,r)) 

x y2h{xv+\i) r -d-l+2h(x +\£) d £ dj ,_ 

Lemma 12.21 allows us to split the integral into $ = $i + (<3? — <l>i), where 
&i(H, A, xo, x, t) is equal to 

( e it\- H i>{x£,r) _ A X 2h(x +X0 r -d-l+2h(x +XO 
^(O.A- 1 ) ^ ' 

Then, 



/ 



|*i(fl-,A,x 0> x,t)| < A 2m / 



xl+ 



it\~ H ip(x,(,,r) _ 



r -d-l+2m ^ 



We notice that 



it\- H ip(x,.) 



< 



2 VOe,.)^) < 2 IV>(x,.)| 



and recall that ip(x,.) belongs to L 1 ^^ x R + ,r _d_1+2m <ir<i£) by Lemma 
Oso that 

lim 4>i(il, A, x ,x,t) = 0. 

A — *0~^~ 

The second term $2 := — <3?i is given by 

(28) -it [ X- H+2h ^ + ^^(x,^,r) r- d - 1+2h ^ +x ^ d^dr . 

JR d x(0,A- 1 ) 

We will now analyse on this expression. 

First, let us note that in the case where h is constant, then the increments 
of X are equal to the increments of X by stationarity and hence $2 = 0. 
That concludes the proof for the fractional microball model, i.e. the first 
part of Theorem 13.11 

Now we deal with the multifractional model (case h non constant). 

a. The smooth case: We assume that (|2*U)) holds and we prove that the 
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critical fdd-lass index of X at xq is > min(l,/3 + 2h(xo)). Without loss of 
generality, we can assume that j5 + 2h(xo) < 1. For H < j3 + 2h(xo), we 
will establish that $2(^7 A, xq, x, t) tends to when A — > + . Recall that by 
equation (|28|) &2(H, A, xo, x, t) is equal to 

—it X~ H + 2h M / ^ r j2(h(xo+X£)-h(xo))^^ x ^ ^ r \ r -d-\+2h(x ) _ 

yR d x(0,A- 1 ) 

Since J R d x ( a- 1 ) ^(^j £' r ) r _d_1+2fe ^°) d^dr vanishes, we get 
$ 2 (^,A,x ,x,t) = -it \- H + 2h ^) 

f {Xr) 2(h( X0+ XO-h(x 0) ) _ A ^( X) ^ r ) r -d-l+2ft(x ) ^ 
^(O.A- 1 ) V ' 

By the /3-Lipschitz assumption on h, for A small enough compared to x, 
using (|2i|) and since = ip 2 , we get 



( Ar )2A(a. o +A0-2fc(*o) _ X j^(x,£,r)| r -*-l+2fc(*o) < C(x)A /3 . 

x(0,A-!) ' 

This implies the lower bound for the fdd-lass index. 

b. The singular case: We now assume that h is /3-Lipschitz on the sphere 
and only deal with the case xq = 0. 

b-i) Assume that {£; h(£) = m} has positive measure. For H = 2m, using 

Lemma E21 again, we get 

(29) 

lim <£> 2 (2m, X,0,x,t) = -it / l h ^ )=m ip(x, £, r ) r -^ 1 + 2m d£dr . 

Hence the finite dimensional distributions of A~ 2m AoAL (A.) converge to the 
finite dimensional distributions of the deterministic field Z m given by 

Z m (x) = - [ l h{0=m Hx, £, r )r~ d - 1+2m <%dr , x G R d . 

iR d xR+ 

It remains to show that Z m is not zero. This follows from the next lemma, 
since {^; h(£) ^ m} contains a ball by continuity of h. 

Lemma 6.1. Let m € (0,1/2). For all Borel sets E C M. d with positive 
measure, if 

( l E {i)^{.,i,r)r- d - l+2m didr = Q 

■/R d xR+ 

in R , t/ien the set E c does not contain any open ball. 
Proof. Note that for all x G M d , 

/ V(^^>~ d ~ 1+2 "X* = (2m -d)" 1 (|x-£|" d+2m - |Cr d+2m ) 
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and hence, for all Borel sets E C 



W(x): = [ l E (OHx,tr)r- d - 1+2m d^dr 
JR d xR+ 

= f l E (0(2m - d)- 1 (\x - £\- d+2m - \£\- d+2m ) d£ . 

jR d V ' 

Let us suppose that we can find an open nonempty ball B C E c . Then Ig 
is smooth on B. The Laplacian of Ig can easily be computed 



JR d 



■d+2m-2 



di 



and proved to be negative on B. Thus, does not vanish on B. This 
completes the proof of Lemma 16.11 □ 

b-ii) Assume that {£; h(£) = m} is of measure 0. we will establish that 
the fdd-lass index of X at is still equal to 2m. For H < 2m, from l|29|l. 
&2(H, A, 0, x, t) tends to 0. On the other hand, for H G (2m, 1) and H < 
2max/i, we have to prove that there exists at least one x S M. d such that 
<&2{H, A, 0, x, t) does not tend to 0. First, let us remark that, for all x € M. d , 
§2{H, A, 0, x, t) may be written as 

-it ( f \- H+2h ^^{x,i,r) r ~ d - 1+2h ^ d^dr + O (X 1 ' 11 )] . 

\JR d xR+ ) 

Thus it is sufficient to consider 

$ 2 (i?,A,x): = f \- H+2h ^^{x,^,r) r~ d - 1+2h ^ d^dr 

JR d xR+ 

(2h(0-d)- l \- H+2h ^f(x,Odt, 



with 

f(x,0 = \x-S\- d+2h ® -\£\- d+2h ®. 
The function /(.,£) is smooth on M. d \ {£} and its Laplacian given by 

Af( x , o = (2h(o - d)(h(o -i)\ x - er d+2h («)- 2 , 

is negative. Then, approximating the Laplacian by the second order in- 
crements, one can find C € (0, ^) such that, whenever \x — £| ^ and 
5<C\x-i\, 

(30) f(x + 8e j ,0 + f(x-6e j ,Z)-2f(x,0<^Af(x,S). 

l<j<d 

Let us take H such that 2m < H < 2 max h and H < 1 and note that the 
sets {£; H < 2/i(£)} and {£,; H > 2/i(£)} have positive measure. Then, by 
continuity of h, there exists a nonempty open ball B C {£; H > 2/i(£)}. For 
every x € R d we introduce 

* A (x) := / 1^(0^(2^(0 " rf)" 1 ^-^®/^^) 
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such that 

* 2 (H,\,x) = * A (x) + (* 2 (-ff, A,x) - $a(x)) . 

By Lebesgue's Theorem, the second term tends to with A. Suppose that 
&\(x) tends to with A for every x in R d . Then for all x G M. d and all 5 G R, 

(31) A< 2) * A (a;) := £ ($ A (x + <5e,) + $ A (x - <%) - 2$ A (») — » . 

l<i<d 

For x & B, and (5 > such that B far, ^) C 5, according to (j30J), 

Af < 5 1 [ (h(o - i)i 2H o< H \x - a\- d+2h ^- 2 d£ < 0. 

Then (|31[) implies that {£; 2/i(£) < i?} has measure zero, which contradicts 
the assumption H > 2m. 

The proof of the first part of Theorem 13.41 is now complete. 



6.2. Fdd-lass properties for the X-ray transforms. Finally we con- 
sider the fdd-lass property at point yo for the X-ray transform. As previ- 
ously, we restrict the computation to the one-dimensional distribution and 
denote V a X = V a X — E (V a X) the centered version of V a X. For any 
y G< a >-*-, tel and H G (0, 1), we write 



Eexp [it\- H (V a X(y + Xy) - V a X(y ))J = exp A, yo, y, t) 
where $>(H, A, yo, y, t) is given by 

f ( e itt-*G,(Ai,,C-i»,r) _ ! _ itX~ H G p {Xy^ - y ,rj) du h {i,r). 

JK d xR+ v 7 

By the same change of variable as in the covariance part of the proof, 
$>(H, A, yo, V, t) is equal to 

A" 1 U^G^iv^r) _ 1 _ it \i- H G p(p+x . ) (y, 7 ,r)) 

<o> J -xRx(0,A- 1 ) v 7 

x (\ r fh(^+pa+yo) r -d-i dldpdr _ 

Let us remark that 

\X-H r A ^o + n 

A ( - T p(p+A.) ► u 

and so 

r i^-«^ + , ) _ 1 _ itt i-H G! ^ A) j „ -^A- 2H+1 G 2 p 2 , 

\ / A — 2j 

where /(A) ~ g (A) if ^M- — ► 1. Consequently we can argue along the 

A^0+ 9[ ' A^0+ 

same lines as in the covariance part of the proof to get, for H = m(a, yo) + 
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1/2, 

*(H,\,y ,y,t) - \t 2 ^l Kyo+pa)=m{a ,y o) p{p) 2d p) r m(a ' W))+1/2 (y,y), 

which concludes the proof. 
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